Complete structure dependent analysis of the decay P —> l + l 
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Abstract 

We use the Mellin-Barnes representation in order to improve the theoretical estimate of mass 
corrections to the width of light pseudoscalar meson decays into a lepton pair, P — > l + l~ . The 
full resummation of the terms (M 2 /A 2 ) n , (m 2 /M 2 ) n and (m 2 /A 2 ) n to the decay amplitude is 
performed, where m is the lepton mass, M is the meson mass and A ~ m p is the characteristic 
scale of the P — ► 7*7* form factor. The total effect of mass corrections is quite important for n{n') 
decays. We also comment on the estimation of the hadronic light-by-light scattering contribution 
to the muon anomalous magnetic moment in the chiral perturbation theory. 

PACS numbers: 
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I. INTRODUCTION 



The theoretical study of the 7r° and r? (77') mesons decaying into lepton pairs and the 
comparison with the experimental rates offers an important low-energy test of the standard 
model. The situation with these decays became more pressing after recent KTeV E799- 
II experiment at Fermilab [lj] in which the pion decay into an electron-positron pair was 
measured with high accuracy using the K L — ► 3n process as a source of tagged neutral pions 

^no-fad (*r° - e + e-) = (7.49 ± 0.29 ± 0.25) • KT 8 . (1) 

The theory for the decay of pseudoscalar mesons to a lepton pair is known for decades 
2, 3, 4], 5]. The main limitation for realistic prediction of these processes comes from the large 
distance contributions of the strong sector of the standard model where the perturbative 
QCD theory does not work. However, it was shown in 6j that theoretical uncertainty can 



be significantly reduced by using CELLO and CLEO data [7, s| and QCD constraints on 
the transition form factors P — ► 7*7*. As a result, the standard model prediction gives 6| 

^Thcor ^0 ^ e+e -) = ^ 2 ± . 1Q - 8j ( 2 ) 

which is 3.3cx below the KTeV result (CQ). It is extremely important to trace possible sources of 
the discrepancy between the experiment and theory. There are a number of possibilities: (i) 
problems with (statistic) experiment procession, (ii) inclusion of QED radiation corrections 
by KTeV is incomplete, (iii) unaccounted mass corrections are important, and (iv) effects 
of newphysics. At the moment, the last possibility was reinvestigated. 



in Q, the contribution of QED radiative corrections to the zr" — < + < decay, which must 
be taken into account when comparing the theoretical prediction (12]) with the experimental 

n n 

data, was revised. In comparison with earlier studies [lOj the main progress made in [9|] 
consists in detailed analysis of the 7*7* — > e + e~ subprocess and revealing the dynamics of 
long and short distances. Occasionally, the final result agrees well with earlier prediction 



based on calculations lOj and, thus, the KTeV analysis of radiative corrections is confirmed. 



There are quite few attempts in the literature to explain the excess of the experimental 



data on the 7r° — ► e + e decay over the standarc 



physics beyond the Standard Model. In Ref. 11|, it was shown that this excess could be 



model predictions as a manifestation of 



explained within the currently popular model of light dark matter involving a low mass (~ 10 



2 



MeV) vector bosons which presumably couples to the axial-vector currents of quarks 



and leptons. Another possibility was proposed in Ref. [12|, [13| in interpreting the same 
experimental effect as the contribution of the light CP-odd Higgs boson appearing in the 
next-to-minimal supersymmetric Standard Model. However, there are other extensions of 
the Standard Model deserving studies in this context. In particular, supersymmetric models 
with R-parity violation and leptoquark models suggest contributions to this process which 
under certain circumstances may be significant. Specific examples of such contributions are 
related to the exchange by the t-squark and the SU 2 l singlet leptoquarks with the couplings 
to u,d-quarks and electrons not stringently constrained from other known processes 14], Il5 |. 

In the present paper, we focus on the mass corrections to the width of light pseudoscalar 
meson decays into a lepton pair, P — > . We show that these corrections are under 
theoretical control and do not help resolving the problem of disagreement of the KTeV data 
on the 7T° — > e + e~ decay with the standard model prediction. However, the mass corrections 
are quite important for realistic predictions for 77 (7/) decays to a lepton pair. 




FIG. 1: Triangle diagram for the P — > l + l process with the pseudoscalar meson form factor 
P — > 7*7* i n the vertex. 



In the lowest order of QED perturbation theory, the decay of the neutral meson, P(q) — > 
l~(p~) + Q 2 = M 2 , p 2 ± = m 2 , (M meson mass, m lepton mass) is described by 

the one-loop Feynman amplitude (Fig. [I]) corresponding to the conversion of the neutral 
meson through two virtual photons into a lepton pair. The normalized branching ratio is 
given by Q, H, Q, 3 



Ro(P -> i + r 



B (P -> l+l' 
B{P^ll) 



(3) 



where (3 (q 2 ) = a/1 — 4m 2 /q 2 and the reduced amplitude is 
d^k (qkf - q 2 k 2 



m 2 (k 2 + ie) [(q — k) 2 + ie] — k) 2 — m 2 + ie] 



p 7 * 7 



(~k 2 , -(q - k) 2 ), 



(4) 



with the transition form factor Fp 7 * 7 * (— k 2 , —q 2 ) being normalized as Fp 7 * 7 * (0, 0) = 1. It is 
convenient to introduce the mass ratio parameters: 

x = (m/A) 2 , u = (m/M) 2 , z = x/u = (M/A) 2 , (5) 

where A ~ M p is a parameter characterizing the transition form factor. In the physically 
interesting cases the parameter x is very small: (m e /A) 2 < 10~ 7 , (m M /A) 2 < 10~ 2 . At the 
same time, we consider the amplitude in the leading order in the fine coupling constant 
a ~ 10~ 2 . Thus, in the following it is reasonable to neglect the x power dependence of the 
amplitude keeping only the dependence on the u and z variables 1 . The aim of this work 
is to improve the previous calculations of the amplitude A(M 2 ) of the P — > decay by 
taking into account all order mass corrections. 



II. RESUMMATION OF POWER CORRECTIONS TO THE AMPLITUDE 

We evaluate the amplitude A (q 2 ) following the way used in ay , • Let us transform the 
integral in (J3J) to the Euclidean metric ho — > ik^. The corresponding integral is convergent 
due to decreasing of Fp 7 * 7 * (k 2 , (q — k) 2 ) in the Euclidean region. It is convenient to introduce 
the modified Mellin-Barnes transformation for the meson form factor 

F Prr (k 2 ,(q-k) 2 ) = -^ [ dz l dz 2 ^(z 1 ,z 2 )r(z 1 )T(z 2 )(^) '(-J*) \ (6) 

[2ixi) Ja+iB? V k J \{k- q) J 

where A is the characteristic scale for the form factor, the vector a = (o"i,o"2) £ ^ 2 , and 

$ (z\, z 2 ) is the inverse Mellin-Barnes transform of the form factor 



-| /*oo /*oo 



(tl,t 2 



(7) 



which is an analytical function of Z\ and z 2 . Introducing Feynman parameterization in (j4j) we 
convert the k— loop integral into the integrals in Feynman parameters which can be expressed 
in terms of T-functions. Then we obtain the following Mellin-Barnes representation for the 
amplitude 16] 



dz\dz 2 dz~i 



-x 



-Zl-Z2„,-Z3 



u Z3 $(zi,2; 2 ) 



2 + —){z 1 + z 2 + z 3 ) 
2u 



CT-HR3 {2nif 

r (-z 3 ) r (i + zi + z 3 ) r (1 + z 2 + z 3 ) r { Zl + z 2 + z 3 ) r (1 - 2 ( Zl + z 2 + z 3 )) 

ziz 2 T (3 - z\ - z 2 ) 



(8) 



Partially, the dependence of the amplitude A (q 2 ) on the x parameter was taken into account in jl 61 ] . 



with a in the 3-dimensional region in space of real parts of Z{ chosen so that the integration 
path a + iM. 3 does not intersect the T-function singularities. Finally, we are able to expand 
the integral over the x, u and z mass ratios by closing the contours of integration in the 
appropriate manner and summing up the obtained series. Neglecting the dependence of the 
amplitude A (q 2 ) on the powers of small x parameter (but keeping \nx) we arrive at the 
following representation: 



A(M 2 ) 



1 



(3{u) 

7T 



2(5 {u 



'1 7T 2 

-In 2 (y (u)) + — 
In (y (u)) 



3, , . 5 3 

+ - In (x) h - 

2 V ; 4 2 

+ A Z , 



dt In 



Li 2 (-y («)) 



t 

A 2 



F P %(t,t) 



(9) 



where 



y(i 2 ) 



P(q 2 



l + (3(q 2 ) 



and the correction to A (M 2 ) is expressed as 
A g = (lux + J ^ (1 - s) 2 [1 - F Prr (sz, 0)] 



1 F W 
2 p 7 * 7 * 



(0,0) 



+ 



1 ds 



Fp 7 » 7 * (t, t) - F Prr (0, t) - tF p ^ r (0, t) 



dy 



(1-3/) 2/ 



1 - 2y + y 3 + ~yFp 7 * 



-s*,0) 



1 - -« 
2 y 



p 7 * 7 * 



-syz, 0) 



1 

+ - 



*d8 
S 



dt \n(t) 2F&2> (t, 0) - Fgw (t - sz, 0) - Fg» it 



(i,o) 



,(1,0) 



(10) 



p 7 * 7 



p 7 * 7 * 



p 7 * 7 * 



— — / ds (1 — s) 2 / dvcosv [ dt -— Im Fp'*K (t — + 2i (stz) 1 ^ 2 cosv, t 

n Jo Jo Jo (stz) ' v 

In the above expressions Fp'*^*(s,i) denotes the derivatives of an order of a and (3 in the 
corresponding arguments of the form factor. 

The first and second lines in are the structure independent parts of the amplitude. 
These expressions agree with obtained earlier results. In particular, they a gree with the 
results of the dispersion approach of 17] and the chiral perturbation theory [rsl. The 
derivation of the amplitude by using the dispersion approach tacitly assumes that the imag- 
inary part of the off-shell amplitude A(q 2 ) is the second line of (jUJ) with M 2 substituted 



by q 2 . Thus, it takes into account only the structure independent part (dependence on the 
parameter u) and is insensitive to the details of the transition form factor. The same result 
appears in the framework of the leading order of the chiral perturbation theory [if], fl9 [- 
because at this order one also does not take into account the form factor. The integral in 
the second line of (JHD was estimated in |6] by using CELLO and CLEO data on the meson 
transition form factors and constraints following from operator product expansion (OPE) 
in QCD. Our new result concerns the structure dependent part of the amplitude A z as the 
combination of the integrals of the meson transition form factor over both the spacelike and 
timelike regions. 

If we retain only the terms linear in z we get 



-\ { (-2 In x + pOft. (0, 0) + 3 j~ dtF P % 



(t,t) 



(11) 



dtHt) Fpjtlf (t, 0) + FpJ^lf* (t, 0) 



(2,0) 



+ 



dt 

~P L 



F Prr (t, t) - F Prr (0, t) - tF P l ;%, (0, t) 



In order to estimate the correction A z to the amplitude related to the meson mass, we 
consider the simplest parameterization of the transition form factor given by the naive vector 
meson dominance model 



7 VMD 

P7*7* 



My 




(M v + s) (M v +ty 

and the correction is 

2 



(12) 




2z - 1 



2z 2 
3 
4? 



In (1 - z) [(1 - z) In (1 - z) + 6z - 4] 



Li 2 (z) + 



z 3 
6 + 2 



(13) 



-arctg 



V4l 



2 - Az + z 2 



+ 



A- z 



(8 — 6z + z 2 } arctg 



V4:Z - z 2 (2 - z) 



2-Az + z 2 



which in the linear approximation (z << 1) reduces to 



A^- 



VMD 



-21nx + - 
3 



(14) 



Taking My = 770 MeV the resulting branchings are given in the Table and compared 
with existing experimental data. The so-called unitary bound appears if in (j3J) only the 



TABLE I: Values of the branchings B {P — > l + l ) obtained in our approach and compared with 
the available experimental results. 





Unitary bound 


CLEO bound 


CLEO±OPE 


This work 


Experiment 


R (vr° -> e+e") x 10 8 


> 4.69 


> 5.85 ± 0.03 


6.23 ±0.12 


6.26 


7.49 ± 0.38 [1] 


Ro(v^ M + M~) x 10 6 


> 4.36 


< 6.23 ±0.12 


5.12 ±0.27 


4.64 


5.8 ±0.8 [20, 21] 


i?o (r? -> e + e~) x 10 9 


> 1.78 


> 4.33 ± 0.02 


4.60 ± 0.09 


5.24 


< 2.7 • 10 4 [22] 


i?o (»/ -> M + aO x 10 7 


> 1.35 


< 1.44 ± 0.01 


1.364 ±0.010 


1.30 




i?o (t/ -> e+e") x 10 10 


> 0.36 


> 1.121 ±0.004 


1.178 ±0.014 


1.86 





imaginary part of the amplitude which is model independent is taken into account. The 
CLEO bound corresponds to the estimate of the real part of the amplitude basing on the 
CELLO and CLEO data on the meson transition from factors. The fourth column of the 
Table contains the predictions where in addition the constraint from OPE QCD on the 
transition form factor for arbitrary photon virtualities is taken into account [6]. 

As expected, a visible change occurs only for r\ (7/) meson decays. It is interesting that 
for 77 decay to muons the mass correction shifts the theoretical prediction in the direction to 
the unitary bound and thus opposite to the experimental result 2lj. This is because the real 



part of the amplitude for this process taken at the physical point q 2 = M 2 remains negative 
and a positive shift due to mass correction reduces the absolute value of the real part of the 
amplitude |Re^4. (M^)|. In this situation, new measurements of 77 — > would be very 

desirable. 

For 7]' decays there appear new thresholds in addition to the two-photon one. In general, 
this violates the so-called unitary bound because the correction A z gains an imaginary part 
at z > 1 



Aim .4 



, 



- ] In 

z 



1 + P 
1-/3 



e(*-i). 



(15) 



As it is seen from the table, this happens for the 7/ — > fifi channel. Nevertheless, it turns out 
that the predictions for this channel are quite accurate. We checked that more sophisticated 
models for the transition form factor (the generalized VMD 23( or the effect of meson mixing 



241 ] ) do not significantly change the results of the Table. 
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III. HADRONIC LIGHT-BY-LIGHT CONTRIBUTION TO MUON G-2 IN CHI- 
RAL PERTURBATION THEORY 



In 25|, the hadronic light-by- light contribution to muon g — 2, which is enhanced by 
large logarithms (LL) and a factor of iV c (~ O (a 3 N c mz) LLj , was estimated in the chiral 
perturbation theory as 

(16) 



LbL,hadr LbL,hadr . LbL,hadr . ^ 

a fi — tt At,LL "+" "/n.pionloop ° ' 



where 



LbL,hadr 



3 /n\ ! /m M ^ViV c V/l ln2 /m| 



16 Vvr 



3tt 



X(A) 



/ (y) = l -y 2 In y - i (2y + 13) + ~ (2 + y) Vy (4 - y) arccos , 

the constant C absorbs subdominant structure dependent contributions. The low energy 
constant \ (A) is related to the amplitude for pion decay into an electron-positron pair by 2 

(17) 



X (A) = 2(31n( I f ) -2.4(0) -7 



-9-6jHctt]n (^)i^ (M)- 



Taking A = 1 GeV we obtain for the logarithmic enhanced term 



LbL,hadr 



(3.4 ±2.0) X 10 



-10 



where uncertainty arises from determination of the integral in (|T 

LbL.hadr 



was a, 



(18) 



The previous estimate 



„ ii -- (5-7ig;°) x 10 — 10 125] . Better accuracy in ([TBI) is due to higher quality of 
determination of x (lGeV) = —17.35 ± 1.20 by using CELLO and CLEO data for the 
transition from factor Fp T * 7 » (t, t) instead of using the decay r\ — *> as was done in 25] 

with the result x (lGeV) 

being of the order 0(a 3 N®) arises from the three-loop 



The second term in (fTST) a 



■14it 

LbL,hadr 
^/i,pionloop 



graphs with a charged pion loop. It was computed in [27j, |28j a^otioop = — 4.46 x 10 



1-10 



2 The logarithmically enhanced hadronic light-by-light contributions to a M are renormalization scheme in- 
dependent. However, the values x(^), / (y) and the constant (1/2) appearing in the renormalization 
group equation depend on the choice of a scheme. Thus, the constant \ (A) defined in [25| is related to 
the corresponding constant \* (tti p ) in 26[ by relation x (A) = —4 [%* (m p ) — 3ln (Jj^-j + l] • 
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Thus, there are large cancellations between In 2 and In terms in cfrj^^* as well as between 



LbL,hadr -i LbL,hadr i r> i i x 

and a^ ipionloop with final result 



a LbLMdr = /_ 1 Q6 ± 2 q + 31C \ x 10 -10 (lg) 

The conclusion is that in the case of P — ► 11 decays the large logarithmic terms are 
dominant and the structure dependent part of the amplitude is a small correction. That is 
why the predictions for the P —>■ 11 decay branchings are rather stable. This is not the case for 
the hadronic light-by-light contribution to muon g — 2. Due to cancellations, the nonleading 
structure dependent terms become crucially important and the chiral perturbation theory 
does not provide realistic estimate. It is necessary to use specific models to provide correct 



predictions for a Lt>L,hadr ^ gee f Qr rev j ew 29J, |3pJ, |3l|, l32| ) . 



IV. CONCLUSIONS 

Our main conclusion is that the inclusion of radiative and mass corrections is unable to 
reduce the discrepancy between the standard model prediction for the 7r° — > e + e~ decay rate 
and experimental result ([T]). Mass corrections are important to get a reliable prediction 
for branchings of 7/-meson to a lepton pair. Moreover, they are small enough to be sensitive 
to details of the transition form factor. For 7/ decay modes not only mass corrections but also 
threshold behavior of the transition form factor are important. Thus, theoretical error for 
T)' decays is less under control comparing with ir° and 77 meson decays. Further independent 



experiments at KLOE [33|, WASAatCOSY 34|], BES III [35j and other facilities will be 
crucial for resolution of the problem. 
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